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1. Introduction 

Neutrino oscillations have been firmly established in the last twelve years or so by 
a beautiful series of experiments with neutrinos from the sun [IH^, the Earth's 
atmosphere [8l[9], nuclear reactors [lOlIlT], and accelerators [12l|13]. All these data 
can be described within a three-flavour neutrino oscillation framework, characterised 
by two mass-squared differences {Aml^, Aml^), three mixing angles (6*12, 6^13, 6'23), and 
one complex phase (5); see section [2?T] for definitions. We know that two out of the 
three mixing angles are large [H], 

sin^ 9,, = 0.3181°:°}^ , sin^ 9,, = O.bOt'o'ol ■ (1) 

The mass-squared differences are determined relatively accurately from the spectral data 
in the KamLAND [TT] and MINOS [H] experiments, respectively |14j . 

Ami, = X 10-^ eV^ , \Aml,\ = 2.40+°:}2 x lO'^ eV^ . (2) 

The parameters in eqs. [1] and |2] are responsible for the dominating oscillation modes 
observed in the experiments mentioned above. 

The topic of this review is the third mixing angle, ^13, whose value is not known at 
present, and is constrained to be small compared to the other two angles [H] (updated 
as of May 2010), 

sin2 0i3 < 0.031 (0.047) 
sin^ 2^13 < 0.12 (0.18) 90% (3(t) CL . (3) 

^13 < 10.1° (12.5°) 

An important contribution to the bound on 9i3 comes from the non-observation of 
disappearance of reactor electron anti-neutrinos at the scale of Am'^, at the CHOOZ [15] 
and Palo Verde [16] experiments, while the final bound is obtained from the combination 
of global neutrino oscillation data, see e.g. [lll[T71[18]. The present information on 6^13 
is reviewed in some detail in section [3l including also a discussion of possible hints for 
a non-zero value [T71IT9] . 

Maybe besides the determination of the absolute neutrino mass and the search for 
lepton number violation in neutrino-less double beta decay, the determination of 6*13 is 
one of the next primary goals in neutrino physics. Its value is of great phenomenological 
as well as theoretical interest. On the phenomenological side, the possibility of CP 
violation in neutrino oscillations, which is a genuine three-flavour effect, depends on a 
non-zero value of ^13. Any realistic possibility to determine the type of the neutrino mass 
hierarchy (i.e., the sign of Am|^) relies on a not-too-small 9i3. Therefore, the results 
on ^^13 from the upcoming generation of experiments will be of crucial importance for 
a possible subsequent high-precision neutrino oscillation facility. We briefly comment 
on the issues of CP violation and mass hierarchy determination in sections 12.21 and 
12. 3[ respectively. Implications of the value of ^13 for neutrino mass models and flavour 
symmetries are briefly mentioned in section 12.41 

There are several neutrino oscillation experiments currently under construction, 
which are expected to start data taking soon. These are the reactor neutrino experiments 
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Daya Bay ^0], Double Chooz '^T\, RENO [22] and the accelerator experiments 
NOi/A [23] and T2K [23]. The primary goal for all of these experiments is the discovery 
of the yet unknown mixing angle ^13. Section H] is devoted to a description of these 
five experiments, and in section [5] we discuss the related phenomenology, highlighting 
the different nature of the experiments (i>e disappearance in reactors versus z/^ — )■ or 
i^^J. appearance in accelerators) and provide sensitivity estimates. 

The results of these experiments will be essential for the planning towards a 
possible next generation of long-baseline neutrino experiments able to address leptonic 
CP violation and the neutrino mass hierarchy. This could be an upgraded super 
beam experiment with a huge detector, or experiments using a new source for an 
intense neutrino beam based on decaying particles in a storage ring, such as a neutrino 
factory [25] (using muons) or a beta beam [261428] (using radioactive ions). These options 
are under intense study, see e.g., [291133] and any decisions will be crucially influenced 
by the results of the upcoming experiments discussed in this review. We briefly discuss 
such future high precision facilities in section El 

We focus in this paper on upcoming experiments which are under construction or 
(from current perspective) are very likely to be funded. Without going into any details, 
let us just mention here also other more speculative ideas towards a 6'i3 measurement. 
These include a high- intensity tritium source immersed into a large TPC [311 135] or 
the use of Mossbauer neutrinos [36H3n] . Such approaches would lead to oscillations 
at baselines of order 10 m, thanks to the low energy of neutrinos emitted in tritium 
decay (Q = 18.6 keV). Unfortunately serious practical as well as principal problems 
seem to make such experiments unlikely from present perspective. The observation of 
neutrinos from a supernova explosion might also allow some conclusions on 6'i3, in lucky 
circumstances even down to values of order 10~^, see for example [1T1 - H3] . Apart from 
the explosion of a nearby supernova this requires one or even more suitable detectors 
at the correct location on earth with respect to the direction of the supernova. ^13 in 
the context of leptonic unitarity triangles has been discussed in [H], implications for 
neutrino-less double beta decay have been considered in [35], discussions about short 
and long baselines electron (anti)neutrino disappearance signals have been published 
in [36]. Implications of ^13 for ultra-high energy neutrinos searched for in neutrino 
telescopes have been discussed e.g., in [3ZII3H], and ^13 effects for neutrinos from WIMP 
annihilations in the centre of the sun have been considered for example in [^[50]. 

Throughout this paper we restrict ourselves to the simplest unitary three-flavour 
framework, and ignore the possibility of additional sterile neutrinos as well as other 
new physics such as non-standard neutrino interactions. Let us note that effects of 
sterile neutrinos are strongly constrained [5T1[52] and one does not expect that they 
have an impact at the sensitivity level of the next generation of experiments discussed 
here. This might be different for a subsequent generation of high precision experiments. 
While model independent bounds on non-standard neutrino interactions are at the level 
of 10~^ — 10~^ [53], in typical theories beyond the Standard Model one expects them to 
be much smaller [531155], beyond the sensitivities of the experiments considered here. 
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2. General remarks on phenomenology 



2.1. Notations and conventions 

Three-flavour lepton mixing is described by a 3 x 3 unitary matrix U, the so-called 
PMNS mixing matrix [561 157]. ^ basis where the charged lepton mass matrix is 
diagonal we have 

3 



■Mi 



(4) 



Here, Ua (a = e, fi, r) are the (left-handed) so-called "flavour fields" participating in 
charged current (CC) interactions, and z/j {i = 1,2,3) are the neutrino mass eigenfields 
with definite masses mi,m2,m^, respectively. Neutrino oscillations depend on the two 
independent mass-squared differences Amgi, Am^^, with Amfj = mf — m^. After 
absorbing unphysical phases in U by redefining charged lepton fields, one finds that 
U contains three complex phases. Two of these are so-called Majorana phases, which 
appear only in lepton number violating processes (such as for example neutrino-less 
double beta decay) but drop out in neutrino oscillations which depend only on the 
so-called Dirac phase. By convention U is parametrised by the three mixing angles 
^12, ^23, ^13, and the Dirac phase 6 in the following way: 
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(6) 



with 



with the abbreviations = sin^j^, Cj^ = cosOj^, and -DMaj 
a and /3 the Majorana phases. The ranges for the mixing angles and the Dirac phase 
are [58] < 6jk < tt/2 and < 6 < 27r. The evolution of the neutrino flavour state can 
be described by a Schrodinger-like evolution equation with the Hamiltonian 

1 



2E, 



-f/diag(0, Am^i, Aml^)U^ + diag(V, 0, 0) 



(7) 



where V is the effective potential in matter responsible for the MSW matter effect [59]- 
|62] . with V = \/2GFNe where Ng is the electron density along the neutrino path. Eq. [7] 
holds for neutrinos, for anti-neutrinos one has to replace U ^ U* and V — > —V. From 
eqs. [5], [6], [7| one can make the following observations: 

• We see from eq. [H] that \Ue3\ = sin ^13. Therefore, ^13 controls the fraction of the 
electron neutrino contained in the neutrino mass eigenfield corresponding to the 
mass m^. Since Am^i <^ |Am|]^| (see eq. [2]), ms corresponds to the mass state 
separated by a larger gap from the other two mass states mi and m2. Hence, for 
^13 = the electron neutrino does not mix with z/3 and can only participate in 
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oscillations with Aml^. This means z/g decouples from oscillations involving m^, 
i.e., oscillations with Amg^^ ^ ^^32- 

• Let us consider the limit ^^13 and Am^]^/! Am|]^| — )• 0. In this case the 
evolution of solar neutrinos and oscillations in the KamLAND reactor experiment 
are described by an effective two-flavour system of and a combination of and 
i^r, governed by the "solar parameters" 612 and Am2i. Oscillations of atmospheric 
neutrinos and in the K2K and MINOS long-baseline z/^ disappearance experiments 
are described by two-flavour Uf^-Ur oscillations governed by the "atmospheric 
parameters" ^23 and Am"^^. These "atmospheric" oscillations are pure vacuum 
oscillations, as it is easy to see from eq. [7] that the matter potential decouples from 
the relevant 2x2 block of the evolution Hamiltonian in this limit. 

• If ^13 > the electron neutrino participates also in oscillations with Amg^. 
Therefore, one can determine 61^ by looking for transitions at the "atmospheric" 
scale involving the electron neutrino flavour. 



2.2. Leptonic CP violation 

Leptogenesis [63] provides a very attractive mechanism to explain the generation of an 
matter-antimatter asymmetry in the early Universe. An important ingredient for this 
to happen is CP violation in the lepton sector. While in general there is no direct 
connection between the CP violation at a very high scale necessary for leptogenesis 
and CP violation observable in low-energy experiments, an observation of CP violation 
in neutrino oscillations would provide a strong hint in favour of the leptogenesis idea. 
Therefore, there is intense activity towards a high-precision neutrino facility able to 
address CP violation in oscillations [29ti32] . 

CP violation in neutrino oscillations is a consequence of a non-trivial Dirac phase 
S. From the parameterisation eq. [5] it is obvious that 6 becomes unphysical if ^13 is 
zero. Leptonic CP violation will manifest itself in a difference of the vacuum oscillation 
probabilities for neutrinos and anti-neutrinos 



where 



P,^^,, - Pp^^p, = -16 J^p sin sin sin , (8) 



Jo,p = Im([/„if/*2t/^if//32) = ±J , J = S12C12S23C23S13C13 sin (5 (9) 

with -|-(— ) for (anti-)cyclic permutation of the indices e,/i, r. J is the leptonic analogue 
to the Jarlskog-invariant in the quark sector [ST], which is a unique and parameterisation 
independent measure for CP violation. It vanishes if any of the three mixing angles is 
zero. Hence, a non-zero 61^ is a necessary prerequisite for leptonic CP violation. A 
recent review on CP violation in neutrino oscillations can be found here [68j . 



Ois-' phenomenology, present status and prospect 



7 



2.3. The neutrino mass hierarchy 

An important question in neutrino physics is the type of the neutrino mass hierarchy, 
which can be normal (NH, Amg^^ > 0) or inverted (IH, Am|^ < 0). This question has 
important consequences for possible neutrino mass models, and the problem of flavour 
in general. Therefore, the determination of the sign of Am|^ is among the main goals for 
a future long-baseline facility, see e.g. [29]. Here we comment briefly on the importance 
of the value of ^13 for this measurement. 

The most promising way to distinguish these two neutrino mass spectra is to search 
for the matter effect in transitions due to Amg^^. The condition for an MSW resonance 
is 

cos2^^i3 = ±^^, (10) 

where +(— ) holds for (anti)neutrinos. For a given sign of Am^^, eq. [TU] can be 
fulfilled either for neutrinos or for anti-neutrinos. Therefore, finding out whether the 
matter resonance due to Am|]^ occurs for neutrinos or anti-neutrinos will provide a 
determination of the sign of /\m\^. As discussed above the occurrence of a matter effect 
in Am|^ transitions is only possible for a non-zero ^13, and hence the possibility to 
determine the neutrino mass hierarchy via the matter effect crucially depends on the 
observability of 6'i3 effects. This can be done in long-baseline experiments [691 [70] or 
with atmospheric neutrinos [7IH7^ . 

In principle changing the sign of Amg^ has also implications for disappearance 
oscillation probabilities in vacuum, independent of the matter effect [751477] . While 
this effect for muon neutrino disappearance is independent of ^13, in the case of 
electron neutrino disappearance also these methods require a large value of ^13. In 
practice, however, such measurements turn out to be extremely difficult, maybe 
unrealistic [781179]. For completeness we also mention the possibility to discriminate 
NH and IH in non-oscillation experiments, e.g. [80l|8T], or from supernova neutrino 
observations, e.g. [H IIS^HSS] . 

2.4. ^13 and models for neutrino mass 

On the theoretical side, the determination of ^13 will provide important information on 
the mechanism of neutrino mass generation and the flavour structure in the lepton sector. 
Naively one may expect that since two mixing angles are large also the third one should 
not be too small. Considering neutrino mass models without any flavour structure, so- 
called anarchical models, one does expect a value of 6*13 close to the present bound [86] . 
If on the contrary experiments would indicate a very tiny value for ^13 one might wish 
to have a symmetry reason as an explanation. 

The intriguing result that mixing in the lepton sector is very different from the 
quark sector might indicate that a special mechanism is at work to produce the peculiar 
values of the lepton mixing angles. For example, rather symmetric patterns for the 
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mixing matrix are the tri-bimaximal [SZ] or the bimaximal |SS] mixing matrices, 



tri-bimax 



/ ^/2/3 1/V3 \ 

y 1/Vq -i/Vs 1/V2 J 



bimax 



/ 1/V2 1/V2 \ 

-1/2 1/2 I/V2 
\ 1/2 -1/2 1/v^ J 



In both cases S23 = 1/2 (maximal 2-3 mixing) and ^13 = 0. While for tri-bimaximal 
mixing = 1/3 (in perfect agreement with data, see eq. [1]), for bimaximal mixing 
= 1/2, which requires significant corrections to be consistent with the experimental 
value. 

The measured values of the mixing angles may be the result of such regular patterns. 
This could indicate a special symmetry among generations. There is a huge literature 
on flavour symmetries, which we cannot properly account for here. Popular symmetries 
are for example a /x-r exchange symmetry [89l[90] (note that both examples given in 
eq. [TT] fulfil the relation = \Uri\) or the permutation symmetry |91] (usually in 
context of tri-bimaximal mixing), see [92] for a recent review and references. In many 
cases flavour symmetries predict ^13 = in the limit where the symmetry is exact, 
as in the tri-bimaximal and bimaximal examples given above. However, one expects 
corrections to the zeroth order symmetric limit due to breaking of the symmetry, which 
then in general will induce a finite value for ^13, see for example [93|[9l]. Furthermore, 
if the symmetry holds at some high scale, renormalisation group running will introduce 
deviations from the symmetric limit, see e.g. [95]. In [96] a model is discussed which 
allows for a quite sizeable value of 6*13 but still preserves the tri-bimaximal predictions 
for 9i2 and 623. 

An alternative approach to flavour comes from Grand Unified Theories (GUTs). A 
priori a GUT by it self makes no statement about flavour (unless it is complemented with 
additional flavour symmetries). However, since quarks and leptons reside in common 
representations of the GUT gauge group, one obtains relations between quark and lepton 
masses and mixing parameters. For example, a specific model based on the SO(IO) gauge 
group has been discussed in [9^, where the predictions for 6*13 based on GUT relations 
between the Yukawa matrices have been worked out explicitly. Predictions of this kind 
can be confronted with future measurements of di^. 

In Ref. [98j a survey of many neutrino mass models has been performed, with a 
particular focus on the predictions for 6^13, see also table 1 of [99]. Here we just want 
to point out the importance of the ^13 measurement for discriminating among models. 
Certainly an improved determination of ^13 (either establishing a finite value or setting 
a more stringent upper bound) will provide important new information regarding the 
problem of flavour. Further discussion and references can be found in the review jlOOj . 



3. Present status and possible hints for 613 > 

The present information on the value of 613 emerges from an interplay of the global 
data on neutrino oscillations, as illustrated in fig. [H see [T^ITSllTOT] for recent global 




analyses. In the following we discuss the phenomenology of the various different data 
sets from the CHOOZ reactor experiment (sec. I3.ip . solar neutrinos and the KamLAND 
reactor experiment (sec. 13. 2p . atmospheric neutrinos (sec. 13. Sp . and the Ve appearance 
data from the MINOS long-baseline experiment (sec. 13. 4p . In section [375] we summarise 
the status of in the global analysis. 

3.1. The hound from the CHOOZ reactor experiment 

An important contribution to the bound comes, of course, from the CHOOZ 
reactor experiment [15]. Experimental aspects and the phenomenology of ^13 reactor 
experiments will be discussed below in sections 14.11 and |5l The CHOOZ experiment 
observed the Ve flux emitted from the two cores of the Chooz nuclear power plant at a 
distance L ^ 1.05 km. The survival probability is given by 

Pee- I- sin^ 2^13 sin^ + Oia") , (12) 

with a = Am2i/Am3]^. Considering typical reactor neutrino energies of A MeV 

one finds that this experiment is sensitive to the ^13 induced oscillations with Amg^, 
whereas terms due to Amgi are suppressed by the small solar mass squared difference 
and can be neglected as long as sin^ 2^13 > ~ 10^'^. The ratio of the measured and 
the predicted neutrino rates in CHOOZ, averaged over the energy spectrum, has been 
obtained as [15] 

R = 1.01 ± 2.8%(stat) ± 2.7%(syst) . (13) 

The parameter constraint resulting from this measurement is shown in fig. [T]in the plane 
of sin^ ^13 and Am|]^. A less constraining result was reported also by the Palo Verde 
reactor experiment [TB]. A meaning full bound on ^13 can only be obtained by combining 
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the CHOOZ result with the determination of |Am|^| from atmospheric and long-baseline 
experiments. From the combined analysis of CHOOZ + atmospheric neutrino data + 
K2K + MINOS disappearance data one finds sin^^ig < 0.027 at 90% CL [H]. 

3.2. Information from solar neutrinos and the KamLAND reactor experiment 

For solar neutrinos as well as the KamLAND reactor experiment it is an excellent 
approximation to consider the limit Am|^ — ?■ oo. Then it follows from the Hamiltonian 
in eq. [7] with the parametrisation for U from eq. |5l that the survival probability of 
electron (anti)neutrinos is given by |102H104] 

Pee ^ c\,Pr, + ^^3 , (14) 

where P'^^ is a two-flavour survival probability depending on 9i2 and ^m\^^ where the 
matter potential V is replaced by c\^V . The s^g term is tiny and will be neglected in 
the following discussion. 

For the KamLAND reactor experiment the matter potential can be neglected and 
P^^ is just the vacuum probability: 

pKamL ^ ^^^4 _ ^^^^2 ^g^^ ^^^^2 . (15) 

This leads to an anti-correlation of sin ^13 and sin^ 612 [105], see also [TOillTOG] . For solar 
neutrinos one obtains simple expressions for the low and high energy part of the solar 
neutrino spectrum, below and above the MSW resonance in the sun at E]^'^ ^ 2 MeV. 
For low energy solar neutrinos one can adopt the vacuum approximation, and averaging 
over the fast oscillations leads to 

psoiar ^ pQg4 (^1 " ^ slu^ 2^i2^ (low energies) . (16) 

The high energy part of the spectrum undergoes the adiabatic MSW conversion inside 
the sun and P^^ ~ sin^ 9 



'12- 



PeT^^"^ ^ cos^ ^13 sin^ 9^2 (high energies) . (17) 

Eq. [16] for low energy solar neutrinos shows a similar anti-correlation between sin^ 9i3 
and sin^ 9i2 as in KamLAND, whereas for the high energy flux subject to the SNO 
CC/NC measurement [5l|6], a positive correlation of sin^6'i3 and sin^ ^12 emerges. As 
discussed e.g. in |104[|105j . this complementarity leads to a non-trivial constraint on ^13. 

Here we present results from the global solar neutrino and KamLAND analysis of 
Ref. [H] (2010 updated arXiv version 3), where a description of the used data can be 
found. For technical details see also |105pi07] . Fig. [2] (left) illustrates the different 
correlations between 9i2 and ^13 from KamLAND/low energy solar and high energy 
solar data discussed above. The right panel shows the Ax^ profiles for two different 
assumptions on the Standard Solar Model. The curves labeled AGSS09 and GS98 
refer to low and high metalicity solar models, respectively |108j . which imply slightly 
different predictions for the *B solar neutrino flux. The "default" assumption in the 




Figure 2. Left: Allowed regions in the (6*12 -6*13) plane at 90% and 99.73% CL (2 dof). 
The black contours correspond to KamLAND combined with the low energy solar 
neutrino experiments Homestake, the Gallium experiments, and Borexino, whereas the 
blue contours correspond to the high energy solar neutrino experiments SuperK and 
SNO. The red shaded region shows the 99.73% CL region for the combined analysis. 
Am2i has been fixed at the combined best fit point. (Figure kindly produced by 
Mariam Tortola.) Right: Ax^ as a function of sin^ 613 for solar and KamLAND data, 
illustrating the impact of the most recent solar neutrino data (mainly the SNO low 
energy threshold analysis [6]) as well as the impact of different solar models (AGSS09 
and GS98 refer to low and high metalicity solar models, respectively |108j ). 



global analysis of [H] is the AGSS09 low metalicity solar model. For a discussion of 
the impact of solar models on the fit see also [18]. Furthermore, the thin black dashed 
curves shows the Ax^ profile for the 2008 solar neutrino analysis from [H], highlighting 
the impact of the most recent solar neutrino data, mainly the SNO low energy threshold 
analysis [6] which leads to an improved determination of the total NC event rate with 
the remarkable precision of 3.1% (stat) and 2.5% (syst) at Icr. 

The 6'i3 fit results [14J from global solar data and KamLAND are (best fit and la 
errors) 

sin^ 013 = 0.022tomt (solar + KamLAND) , (18) 

and the following bound is obtained at 90% {3a) CL: 

sin^ 613 < 0.053 (0.078) (solar + KamLAND) . (19) 

The result from eq. [18] is in excellent agreement with the value obtained by the SNO 
Collaboration [6], sin^ ^13 = 0.0200+[!:JJ?|?^, and it implies a slight preference for a non- 
zero value of ^13. At ^^13 = we have Ax^ = 2.2 (see fig. [2]), which corresponds to a 
hint at about 1.5cr. This hint emerges from a slight mismatch of the best fit points for 
6112 and Aml^ for 6113 = for solar and KamLAND data separately pT ] nM ][T06l[T07] . 
This mismatch can be alleviated to some extent by a finite value of ^13, see e.g. |106j for 
a discussion. From fig. [2] (left) we see that the high energy solar neutrino experiments 
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SuperK and SNO show a slight preference for a non-zero 6*13. Eq. [T7] indicates a 
degeneracy among 613 and 6*12 as visible in the figure. This degeneracy is not complete 
due to additional information from the spectral shape as well as day-night asymmetries 
induced by the matter effect in the earth during night. While the slope along the 
SuperK-|-SNO allowed region is very flat, these data pull the global best fit point 
somewhat towards a non-zero value. For models with higher solar metallicities like 
GS98, a slightly larger best fit point is obtained, sin^^ig = 0.027l^:|J}^ and Ax^ = 3.05 
at ^^13 = 0, while the bound on 613 becomes slightly weaker, see fig. [2] (right). 



3.3. ^13 and atmospheric neutrinos 

Atmospheric neutrino data come from the SuperKamiokande experiment that recently 
released an updated analysis of the combined SK I+II+II data |109] . Atmospheric 
neutrinos provide a probe of several decades in L/E^ thanks to the wide range in neutrino 
energies, from few 100 MeV to several GeV, and baselines of L ~ 15 km for down going 
neutrinos up to L ~ 12 000 km for up-going ones. Furthermore, up-going neutrinos 
have large trajectories through earth matter, opening the possibility to look for matter 
effects. On the other hand, the atmospheric neutrino "beam" is very messy, consisting 
of u^, Ufj^, z/g, z/g fluxes with large uncertainties. Therefore, from the observation of muons 
and electrons from neutrino CC interactions one can access only certain combinations of 
transition and survival probabilities. Also, the neutrino direction for a given CC event 
is not known. 

In water Cerenkov detectors such as SuperKamiokande, electron-like events provide 
sensitivity to three-flavour effects. One can identify three type of effects: di^- 
induced [7T |lll01fll6] . Am2i-induced |117H121j . and interference effects of the previous 
two |119j . see e.g. |122j for a discussion. Effects induced by ^13 are important in the 
multi-GeV energy range, whereas Am^i effects are mainly relevant for sub-GeV energies. 
Defining the excess of e-like events as eg = {Ng/N^ — 1), with (Nj^) being the number 
of e-like events with (without) oscillations, one has 

er'" ^ (r sin^ ^23 - l) {P^D , (20) 
^ (rcos2e23-l)(P|r)- (21) 
Here (P^i) {{P21)) is an effective two-flavour probability governed by Am^i and 
6*13 (Am2i and 6^12), appropriately averaged and including the weighted contributions 
from neutrinos and anti-neutrinos, and r{Ey,Q) = F^/F^ is the ratio of the initial 
unoscillated muon and electron neutrino fluxes, B being the zenith angle. Since for sub- 
GeV energies r ^ 2, Am\^ effects are suppressed for ^^23 ~ 7r/4, however they provide 
a sensitive measure for deviations from maximal 6'23 mixing. Indeed, small deviations 
from 6'23 = 7r/4 found in |122[I123] may be traced back to a slight excess of sub-GeV 
e-like events in the SuperKamiokande data, see, however |109j . where such deviations 
are not found. 

Here we are more interested in ^13 induced effects from eq. [201 relevant in the multi- 
GeV range, where one has r f» 2.6 — 4.5. Hence one expects that these effects could show 
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up also for maximal 623 mixing. Qualitatively, e'f^"^^^ vanishes for ^13 = and increases 
monotonically with ^^13. The effect is most pronounced for zenith angles corresponding 
to neutrino trajectories crossing the earth mantle, or earth mantle and core, where ^13- 
effects can be resonantly enhanced due to matter effects j59 | l6T |[TTmilll] . In the relevant 
zenith angle bins e™"'*' can reach values of the order of 10% (see e.g. fig. 5 of [71]). For 
the normal hierarchy the resonant matter enhancement occurs for neutrinos, whereas 
for the inverted hierarchy it occurs for anti-neutrinos. Since the event numbers in water 
Cerenkov detectors are dominated by neutrinos because of larger cross sections, e™'^!*' is 
larger by a factor of 1.5 — 2 for the normal hierarchy than for the inverted one. 

In addition to Am^i-effects of eq. ED and 6'13-effects of eq. [20] also an 
interference term between the two contributions is present |119] . It is proportional to 
(r sin ^13 sin 2^23), it vanishes in the limit Amg^ = 0, and it depends on the CP-phase S. 
Because of the different dependence on the flux ratio r the interference term may become 
important in cases where the effects governed by eqs. [21] and [20] are suppressed |122j . 

Such three-flavour effects in atmospheric neutrinos can be used to determine the 
neutrino mass hierarchy or the octant of 623 (if different from 7r/4), in huge future 
atmospheric neutrino detectors (e.g., water Cerenkov or magnetised iron), possibly in 
combination with a future long-baseline experiment, see for example [721471] [T231 - I1 26] . 
In the following we discuss possible implications of present data from SuperKamiokande 
(SK) on 613. 

In Ref. [122] a preference for a non-zero 613 value from SK-I data was noted 
(see also [T27] ). Refs. [T71[T0T| [T22] find from atmospheric + CHOOZ + long-baseline 
disappearance data a 0.9a hint for a non-zero value: sin^ ^13 = 0.012 ± 0.013. In the 
atmospheric neutrino analysis in [M] which is based on [105] Am2i effects are neglected, 
and combined with CHOOZ data the best fit occurs for ^13 = (c.f. fig. [3l right). This 
is in agreement with a similar analysis by the SuperKamiokande collaboration |109U128j . 
Also, in the atmospheric neutrino analysis from Ref. [123] (which does include Arn^i 
effects, as Refs. [T7 ]ll22j ) the preference for a non-zero ^13 is much weaker than the one 
from [122], with a Ax^ < 0.2. 

A discussion and comparison of the results of different groups can be found in |129j . 
The possible origin of the hint from atmospheric data has been investigated in some 
detail in [T8lll30j . Ref. [130] concludes that the statistical relevance of the hint for 
non-zero ^13 from atmospheric data depends strongly on the details of the event rate 
calculations and of the analysis (such as e.g., the specific way of how systematic 
uncertainties are treated). Furthermore, the reason for the hint seems be a small excess 
of multi-GeV e-like events in SK-I data (used in [T7 ]I122] ). which however disappears 
after the inclusion of SK-II data. Ref. [18j comes to a similar conclusion from the 
combined SK-I-I-II-I-III data. 

Let us stress, however, that all analyses agree within Ax^ ~ 1 and therefore 
there is no significant disagreement. Certainly these are very subtle effects 
sensitive to the fine details of the data analysis which can be addressed optimally 
only within the experimental collaboration. The recent three-flavour analysis by 
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SuperKamiokande |109] did not find any hint for a non-zero 6*13. Tlie ^13 analysis 
of |109] lias been performed adopting the approximation Am^i = 0, and therefore, a 
direct comparison with the results of the phenomenological groups pT t [T8 |[T22|ll23lll30] 
is still not possible. 

3-4- Ve appearance data at accelerators 

The first result on appearance at a long baseline experiment has been published by 
the K2K collaboration [m]: sin2 2^i3 > 0.3 (90% CL) in a two-flavour approximation 
analysis. Now also the MINOS experiment provides a first glance at a z/^ — )• Ue 
appearance search for ^13 [19], which is actually the main objective for the upcoming 
T2K and NOz/A experiments, see section HI 

In Ref. p2j a search for z/^ — )■ z/g transitions by the MINOS experiment has been 
presented, based on a 3.14 x 10^*^ protons-on-target (pot) exposure in the Fermilab NuMI 
beam (out of the 7 x 10^° pot accumulated so far). 35 events have been observed in 
the far detector with a background of 27 ± 5(stat) ± 2(syst) events predicted by the 
measurements in the near detector. This corresponds to an excess of about 1.5a which 
can be interpreted as a weak hint for z/g appearance due to a non-zero 6^13. 

In the MINOS detector, being optimised for muons, it is rather difficult to identify 
z/g CC events since they lead to an electromagnetic shower not very different from a 
TT° signal. Neutral current (NC) and misidentified z/^ CC events often have a similar 
signature, and hence lead to a background for the z/g appearance search. Indeed, in 
Ref. |132j an analysis of "NC events" has been performed, where "NC events" in fact 
include also z/g CC events due to the similar event topology. Therefore, a possible 
^'/i oscillation signal would contribute to the "NC event" sample of |132j and these 
data can be used to constraint 6*13. 

We have performed a fit to the MINOS z^g appearance and NC data by using the 
GLoBES simulation software |133[ll34j . The predicted z/g appearance spectrum has 
been calibrated by using the information given in |135] . In the fit we include a 7.3% 
uncertainty on the background normalisation (Tab. I of [I9]). For the NC analysis we 
have performed a fit to the observed spectrum by summing the NC events induced 
from the total neutrino fiux with the z/g CC appearance signal due to oscillations. We 
include a 4% error on the predicted NC spectrum and a 3% error on the z/^j CC induced 
background (Tab. II of |132] ). A full three-fiavour fit is performed taking into account 
a 5% uncertainty on the matter density along the neutrino path. 

The z/g — )■ z/^ oscillation probability depends in a non-trivial way on all 6 oscillation 
parameters, in particular on the unknown value of the CP-phase 6 as well as on the 
sign of Amg^ (neutrino mass hierarchy). The dependence has to be taken into account 
when extracting information on ^^13. This will be discussed in more detail in section [5] 

I We note however that the signal excess depends very much on the choice of the position of the cut 
on the neural network variable that separates signal from backgrounds (fig. 2 of reference [H]). Also 
for this reason the publication of the z/g appearance result with the full statistics is quite important. 
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Figure 3. Left: Allowed regions in the (sin^ 6*13 - S) plane at 68% CL (2 dot) 
for MINOS i^e appearance and NC data. Regions are shown separately for normal 
(NH) and inverted (IH) neutrino mass hierarchy. For comparison we show also the 
bound from global data at 90% CL (1 dof). Right: Ax^ projection as a function 
of sin^ 9i3 for MINOS appearance and NC data, assuming NH (solid) and IH 
(dashed), both with respect to the common minimum, which occurs for IH. The 
green solid curve corresponds to the bound from CHOOZ + atmospheric + K2K + 
MINOS (disappearance) data. 



in the context of the upcoming long-basehne accelerator experiments. Fig. [3] shows the 
fit results for MINOS z/g appearance and NC data, illustrating the importance of the 
parameter correlations. The left panel shows the allowed regions in the (sin^ 9i3 — S) 
plane, for both neutrino mass hierarchies separately. We observe the typical S-shaped 
regions coming from the trigonometric dependence on the phase S. In the right panel 
the is marginalised with respect to all parameters except for 613 and sgn{Am'^i), 
where for the solar and atmospheric parameters we imposed Gaussian errors taken from 
eqs. [H [21 without including any other information on 6*13 except from MINOS itself. 

First, we note that the best fit point is always obtained for the inverted hierarchy 
(Am|^ < 0), and in that case in general the constraint on sin^6'i3 is weaker, since for 
IH the matter effect tends to suppress the z/g appearance probability. In fig. [3] (right) 
the Ax^ for normal hierarchy (Amg^^ > 0) is given with respect to the best fit for IH. In 
the global analysis we also marginalise over the two hierarchies, and hence, the actual 
information from MINOS comes from the IH. 

Second, we see from the figure that MINOS z/g appearance data shows a slight 
preference for a non-zero value of ^13, with a best fit point of sin^ 6'i3 = 0.032(0.043) 
for NH (IH) with Ax^ = 1-8 at sin^6'i3 = (about 1.3cr). In contrast, no indication 
for a non-zero ^13 comes from the NC data. Furthermore, one observes that NC gives 
a slightly more constraining upper bound on sin^ ^13 than appearance, while both 
are significantly weaker than the bound from z/^ disappearance data + CHOOZ. Let 
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us mention that the result for the NC analysis strongly depends on the value assumed 
for the systematic uncertainty, whereas the appearance result is more robust with 
respect to systematics, being dominated by statistics. 

MINOS Ue and NC data are not independent, and adding the corresponding x^'s 
would imply a double counting of the same data. Therefore, only z/g appearance data 
without the information from NC data is used in the global analysis of [H] . It is found, 
however, that adding both MINOS data sets (ignoring the double counting problem) 
leads to practically the same result in the global fit, both for the "hint" for ^^13 > as 
well as the global bound, the latter being dominated by other data sets. 

The OPERA experiment [I36] at the CERN to Gran Sasso (CNGS) beam has 
an excellent electron neutrino selection efficiency, but the beam setup is optimised for 
i^fi^i^T transition searches, being not very competitive for electron neutrino appearance 
searches. OPERA sensitivities to 613 have been firstly computed in |137] and then 
in [T38] where a sensitivity of sin2 2^i3 > 0.14 (90% C.L., for Am^^ = lO'^ eV^ ) has 
been estimated in a 5 years neutrino run at the nominal CNGS intensity of 4.5 ■ 10^^ 
pot/year (expected to happen after 2013). Possible upgrades of CNGS have been studied 
in |139[I140] where an off-axis liquid argon detector would detect neutrinos at the first 
oscillation maximum. For a discussion of these experimental possibilities see |141] and 
also section [61 

Note added: After the completion of this review the appearance data from 
MINOS corresponding to 7x 10^° pot has been made public |142j . 54 events are observed 
with an expected background of 49.1±7.0±2.7. The 1.5a excess found initially in pL9j has 
reduced now to about 0.7cr and data are therefore in perfect agreement with background 
expectations for no appearance. The results and figures in this section correspond to 
the initial 3.14 x 10^'^ pot, while final results on 613 in the global analysis from [ll] have 
been updated with the 7 x lO^*' pot data. 

3.5. 6'i3 in the global three-flavour fit 

Let us now summarise the present situation obtained in the global fit of all relevant 
oscillation data, as illustrated in fig. [H In the updated analysi^ of [H] the following 
bounds at 90% (3cr) CL are obtained (c.f. fig. [H right): 



§ The results presented here are based on the arXiv version 3 of [M], updated with the 7 x 10^" pot 
MINOS Ve appearance resuhs from [ 142] . 
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Table 1. Comparison of the best-fit values for sin^ ^13 and the significance of the 
hint for 6*13 > from different global fits to neutrino oscillation data. The numbers 
from [IH] and [H] include 7 x 10^° pot Ve appearance data from MINOS, whereas |101| 
is based on 3.14 x 10^° pot. AGSS09 and GS98 refer to low and high metahcity solar 
models, respectively |108] . 



The "hint" for ^^13 > coming from the different data sets can be quantified by 
considering the Ax^ for ^^13 = 0: 

solar+KamLAND 
CH00Z+atm+K2K+MIN0S 
MINOS z/g appearance 
global data 

In table [U we compare the best-fit values for sin^ ^13 and the significance of the hint 
for 6'i3 > from the global fits to neutrino oscillation data from three different groups. 
All groups find a non-zero best-fit point in the range sin^ ^13 = 0.01 — 0.02. Depending 
on the analysis as well as variations in assumptions about the solar metalicity used as 
input for the solar neutrino flux prediction a "significance" for a non-zero ^13 between 
l.lo" and 2a is found. Part of the differences can be attributed to differences in the 
atmospheric neutrino analyses, as discussed in section 13.31 However, it is the opinion 
of the authors that within the accuracy which one can attribute to global analyses of 
this kind the results are in agreement. While it is premature to draw strong conclusions 
from these results, upcoming experiments will answer very soon the question whether 
6*13 is indeed in the range indicated by present global analyses. This will be the topic of 
most of the remaining part of this review. 

Let us also mention that there are some limitations of the statistical method applied 
to obtain these results. All three groups quoted in tabled] use a Ax^ method to evaluate 
allowed regions as well as upper bounds on 613 assuming standard x^-distributions. For 
example, the 90% CL bound is obtained by the requirement Ax^ = 2.71. The confidence 
levels obtained by this method are only approximate close to the physical boundary of a 
parameter, such as sin^ 9is > in our case of interest. The confidence interval from such 
a Ax"^ method should be interpreted as a two-sided confidence interval, which lacks a 
well defined meaning close to the boundary. Therefore, the results on ^13 quoted in this 
section should be taken with some grain of salt, and in particular the numbers given 
for various confidence levels (for the significance of a hint for a non-zero 613 as well as 
upper bounds on 6*13) have to be considered only as approximate, and should always be 
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Figure 4. Neutrino induced positron spectrum for various values of sin^ 2^13 simulated 
for the Double Chooz far detector |143) . 



understood in terms of the Ax^ values. To convert the Ax^ values into a well defined 
probability would require more elaborate statistical methods than used so far in the 
literature. 

4. Description of upcoming experiments 

We move now to the discussion of the upcoming generation of experiments searching for 
6'i3, reviewing some general principles of reactor experiments (sec. 14. ip and accelerator 
neutrino beam experiments (sec. 14. 2p . In both cases we present the representatives of 
these classes of experiments: Daya Bay [20], Double Chooz [21], RENO |22j for the 
reactor experiments, and NOz/A [23] and T2K [21] for the beams. 

4.I. Reactor experiments 

Reactor experiments see a large signal of z/g events, and search for a small deviation from 
the non-oscillation prediction due to ^13-induced z/g disappearance, see fig. HI This is a 
precision experiment, whose success relies on statistical as well as systematical errors 
below the percent level. 

The detection of nuclear reactor anti-neutrinos is a rather well known topic, since 
neutrinos themselves have been discovered at a nuclear reactor in 1956 |144] . Modern 
experiments are all inspired by the CHOOZ detector [145j . and this technique is nicely 
described in |146j . Here we will very briefiy summarise the basic principles. 

The signal of a reactor anti-neutrino is provided by the reaction 

Ve + p ^ + n , (24) 

a process with a threshold energy of 1.8 MeV where the positron provides a prompt signal 
(boosted by the two 511 keV annihilation gamma rays) and the neutron a delayed signal 
(in liquid scintillator this delay is about 170 /is). The coincidence of the two signals 
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is powerful enough to make it possible to run the experiment with a tolerable level of 
backgrounds. 

Liquid scintillator can be loaded with 0.1% natural gadolinium, an element with 
high thermal neutron capture cross-section. In this way the neutron capture time 
is reduced from ~ 170 /xs to ~ 27 /xs, allowing for the reduction of the uncorrelated 
background. Furthermore, Gd de-excitation after capture releases an 8 MeV 7 cascade 
(it would be ~ 2 MeV in pure hquid scintillator), producing an integrated signal well 
above the natural radioactivity. 

The CHOOZ experiment concluded its data taking with a 2.8% statistical and 
a 2.7% systematical error (c.f. eq. [T3l) . The goals of a follow-up experiment are to 
improve CHOOZ sensitivity by a factor 5 at least. This roughly reflects on a factor five 
improvement both in statistics and in systematics. 

Let's start on the easy part, the statistics. The CHOOZ detector was a 5 ton 
detector exposed to two reactors of 8.6 MW thermal power at a distance of 1.05 km. 
The experiment integrated a total run of 8761.7 h, only 1543.1 of which with the two 
reactors on and 3245.8 h with one of the two reactors on. To gain a factor 25 in the 
number of neutrinos a detector twice as big running 3 years with an improved efficiency 
is needed. The main limiting factor in this direction is the stability of a gadolinium 
doped liquid scintillation detector. Important progress has been made in this field in 
recent years [147] such that a running time of 5 or more years seems feasible. 

More generally the number of events can be estimated by 



where L is the baseline, T is the running time of the experiment, Pth is the thermal 
power of the nuclear power plant, M is the fiducial detector mass, and we have assumed 
100% efficiency. 

More difficult is the reduction of systematics. This task is particularly challenging 
because CHOOZ had the very rare opportunity of directly measuring the backgrounds 
with the nuclear reactors off. For a complete discussion and comparison of systematic 
errors at reactor experiments see [148] . 

The first important action is to introduce a close detector as similar as possible 
to the far detector |14911152j . in order to measure the neutrino interaction rate before 
oscillations. In this way uncertainties on the neutrino rates (around 2%) almost cancel 
out. This approach has some intrinsic limitations. The neutrino rates can't be the same 
in the two detectors (even in absence of oscillation) because of the different coverage and 
of the different distance from the source. The calibration and live-time of the detectors 
must be kept as similar as possible. In a configuration where more than one reactors 
are present, the close detector doesn't measure the identical fiux of the far one. This is 
because the neutrino fiux of a reactor varies in time according to the core composition 
which differs from reactor to reactor following their fuel refurbishment. 

The backgrounds of reactor neutrino experiments are of two types: uncorrelated 
signals from cosmic rays and natural radioactivity and correlated signals from neutrons 
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Figure 5. Number of disappeared events in Double Chooz for sin^ 2^13 = 0.1, i.e., the 
difference between no-oscillation and oscillation spectra in the far detector, compare 
fig.m Also shown are statistic and expected systematic errors, together with systematic 
errors as big as the former Chooz experiment. The signal is compared to the expected 
background rate. Based on data from [21]. 

generated by cosmic muons. Uncorrelated signals can be separately measured and 
normalised. Their random coincidences can generate background events that can be 
kept at a negligible level with a careful choice of low activity materials and design of 
the detector. 

More problematic are neutrons induced by cosmic muons. A cosmic muon not 
crossing the detector can induce spallation or be captured on the materials outside the 
detector. These can produce neutrons that escape the vetoes and produce correlated 
signals inside the detector that can mimic the anti-neutrino signals. Among the 
cosmogenically produced isotopes there are some long lived, as '^He and ^Li, with 
decay times of 119 ms and 174 ms respectively |153] . that make hardware active vetoes 
impractical. As an example, the close detector of RENO is expected to have a 25% 
dead time for a 0.5 ms veto after any detected muon in the outer veto. This is the main 
reason that forces detectors to run at shallow depth, complicating very much the choice 
of the possible sites. 

Fig. |5] illustrates the size of the disappearance signal in Double Chooz, compared 
to statistical and systematical errors as well as the expected background. Even if 
clearly inspired by the CHOOZ design, next generation reactor neutrino detectors 
have introduced improvements in the design in order to reduce the primary sources 
of systematics: knowledge of the fiducial volume and backgrounds. 

In the following, as an example, we will describe the detector of the Double Chooz 
experiment, shown in fig. [6l The inner detector is made of 10.3 m'^ (8.3 metric tons, 
with a diameter of 2.3 m) of Gd-loaded (0.1%) liquid scintillator in an acrylic vessel. It 
is immersed in a gamma catcher of 22 m^ (3.4 m of diameter) of undoped scintillator, 
aimed to detect the gammas emitted in both the neutron-capture process and positron 
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Pxh [GW] 


L 


111 


"^Det [t] 
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Backgrounds/day 


Daya Bay [20| 


17.4 


1700 


80 


10 


10^ 


0.4 


Double Chooz [2l] 


8.6 


1050 


8.3 


1.5 


10^ 


3.6 


RENO p2J 


16.4 


1400 
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3 
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Table 2. Summary of experimental key parameters of upeoming reactor neutrino 
experiments. We give the thermal reactor power, the approximate distance between 
reactors and far detector, and detector mass, neutrino events per year, and background 
events per day, all for the far detector. RENO backgrounds are the sum of correlated 
backgrounds as computed in |22^ and uncorrelated backgrounds as estimated in [148) . 




annihilation in the target. In this way gammas emitted from signal neutrino events 
in the outer volume of the target are detected, providing a well-defined target volume. 
A third shield of non-scintillating paraffin oil, 5.5 m diameter, separates the active 
target and the gamma catcher from the photo-multipliers, greatly reducing the intrinsic 
radioactivity of the 390 10-inch photo-multipliers, the most radioactive component of 
the detector (this third shield was not installed inside the CHOOZ detector). The outer 
detector volume is steel walled, 6.6 m diameter, filled with scintillator and lined with 70 
8-inch photo-multipliers, equipping the Inner Veto, having the purpose of detecting and 
tracking muons and fast neutrons. An Outer Veto is placed on top of the detector, made 
of strips of plastic scintillator and wavelength-shifting fibres. Its task is to tag muons 
interacting around the detector producing cosmogenic isotopes, some of which produce 
correlated backgrounds in the inner detector, the most dangerous source of backgrounds 
for the experiment. 

Below we describe briefly the three next generation reactor experiments. 
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Figure 7. Configuration of the experimental layout of Double Chooz, RENO, and 
Daya Bay. The dashed curve in the Double Chooz configuration is the far flux iso- 
ratio curve. Updated from 11481. 



Table [2] summarises a few key parameters and fig. [7] illustrates the reactor/detector 
configurations. 

4.1.1. Double Chooz The Double Chooz experiment [2T] is being installed near the 
Chooz two-core (4.27+4.27 GW) nuclear power plant. The far detector, described in 
the previous section, is at 1.05 km from the two reactor cores, in the same site as 
the original CHOOZ experiment, at a depth of about 300 m.w.e. (meters of water 
equivalent). The close detector is designed to be identical to the far detector, it is 
placed 400 m from the two reactors, at a depth of 115 m.w.e., see fig. [71 Let us mention 
that for Amli ~ 2.5 x 10"'^ eV^ the far detector distance of 1 km is somewhat too 
short, shifting the oscillation signal to the lower part of the spectrum, whereas the near 
detector distance of 400 m is somewhat too far, with some effect of oscillations already 
present, see e.g. |152] for far and near detector baseline optimisation studies. 

The number of neutrinos detected in the far detector assuming 3 years running time 
will be ~ 45000 compared to 2700 in the Chooz experiment, reducing the statistical 
error from 2.8% to 0.47%. The goal about systematic errors is to reach a level of 0.6%. 
Without the near detector the systematic error would be about 2.5%. 

4.1.2. Daya Bay The Daya Bay experiment |20] will receive neutrinos from three 
nuclear plants, each consisting of two cores: Daya Bay, Ling Ao I and Ling Ao II 
(scheduled to be commissioned by the end of 2010) located in the south of China, 55 
km to the northeast of Hong Kong. The thermal power of each core is 2.9 GW, hence 
the existing total thermal power is 11.6 GW, and will be 17.4 GW after 2010. 

The basic experimental layout of Daya Bay consists of three underground 
experimental halls, one far and two close, linked by horizontal tunnels under 
construction. The geometry of the reactor cores, the two near detector stations and 
the far detector station is illustrated in fig. [71 It is evident from this configuration that 
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the flux measured by the near detectors will not be identical to the one seen by the far 
detector. Each close detector hall, at a depth of about 100 m.w.e., will host two 20 ton 
gadolinium doped (0.1%) liquid scintillator detectors, while the far hall, 350 m.w.e., will 
host four such detectors. The detectors are designed to be movable so that the close 
detectors and the far ones can be swapped. This swapping is not necessary to reach the 
designed systematics (0.38%), but could be performed to cross check the sensitivity and 
possibly further reduce the systematic errors (down to 0.18%). 

The design of the detectors is very similar to the Double Chooz one, with the 
notable difference that the three stations will be submerged in a water pool of 2.5 m 
depth to shield the detectors from ambient radiation and spallation neutrons, see fig. El 
Above the pool a muon tracking detector made of 4 layers of resistive-plate chambers 
(RPCs) will be installed. 

4-. 1.3. Reno The RENO [22] experiment is located on the site of the Yonggwang nuclear 
power plant in the southwestern part of Korea. The plant consists of six reactors lined 
up in roughly equal distances and spans about 1.3 km. With a total thermal power of 
16.4 GW it is at present the second largest nuclear reactor plant in the world. 

RENO will use two identical detectors, a close detector at about 290 m from the 
reactor array (at a depth of about 110 m.w.e.) and a far detector at 1380 m (at a 
depth of about 450 m.w.e.), see fig. |71 The design of the RENO detectors is very 
similar to Double Chooz, the active target will be 15.5 ton of Gd loaded (0.1%) liquid 
scintillator. The outer veto system is a layer of water, 1.5 m thick, contained in a 30 
cm thick concrete vessel. PMTs are mounted on the inner surface of the veto container 
for detecting Cerenkov light from cosmic muons. Goal of the experiment is to reach a 
< 0.5% systematic error. 
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Figure 9. Muon neutrino interaction rates in MINOS, NOz^A, and T2K assuming no 
oscillations, normalised to one year of operation at the power of 410, 700, and 750 KW, 
respectively. The horizontal dotted lines represent the neutrino energy region where 
the Vfj^ — i^e oscillation probability is bigger than 50% of the maximum {6 = and 
normal hierarchy). The MINOS probability is very similar to NOz^A. 



4-2. Accelerator experiments 

Accelerator experiments look for the appearance of the u^. flavour in an almost pure 
z/^ beam, due to oscillations. The background rate has an intrinsic component given 
by the Uf, contamination in the neutrino beam, ranging from 0.5% to 1%. Detector 
backgrounds can be generated by NC events where a tt° produces a signal misidentified 
as an electron, or CC events where the muon is misidentified as an electron. 

Several experiments already tried in the past to detect —¥ z/g transitions at 
accelerators with short baselines. The most stringent result comes from the NOMAD 
experiment that reached a sensitivity on the probability of these transitions -P(t'^ — )■ 
i^e) < 0.7 • 10~^ (90% CL) |154] . The challenge of the next generation of long baseline 
experiments is to obtain similar sensitivities with detectors 10^ more massive than 
the NOMAD detector |155j . a 3 ton, low density (0.1 g/cm^), very sophisticated 
spectrometer. The pioneering z/^ — )■ z/g long-baseline appearance searches have been 
performed by the K2K |131] and MINOS [19] experiments, see section 13.41 

A rough estimate for the sensitivity can be obtained by comparing the signal to the 
statistical and systematical uncertainty of the background: 



where 5* (B) is the number of signal (background) events, and cxbg is the uncertainty in 
the background. 

It is almost impossible to derive a general rule for the number of events in the far 



S 
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Figure 10. Signal and backgrounds events for T2K, computed for sin 26i3 = 0.1, 
6 = 0, AvTigj^ = 3 • 10^"^ eV^ and normal hierarchy. From [156) . 



detector, since it depends very much on the geometry of the target, the optics of the 
neutrino beam line, and the off-axis angle. The number of signal events S is proportional 
to the detector mass times the power of the primary proton beam. Surprisingly enough 
5* does not depend very much on the primary proton energy and on the baseline L 
assuming the far detector at the first oscillation maximum. This implies a constant 
value of L/E with E the neutrino energy. Since the neutrino cross section generally 
increases as E^ with r ^ 1 — 2, the scaling of the beam flux is largely compensated. 
A comparison of the neutrino interaction rates in T2K, NOi/A and MINOS is reported 
in fig. 1 

The SuperKamiokande water Cerenkov detector used in T2K already demonstrated 
to be capable of keeping detector backgrounds to a level smaller than the intrinsic 
h>e contamination of conventional neutrino beams. Also the total active scintillator 
technology (TASD) chosen for NOi/A allows to keep beam and detector backgrounds 
at comparable level. As an example of signal and background events in an accelerator 
experiment we reproduce signal and background events in T2K as shown in [156] . 

However, the water Cerenkov and TASD techniques allow a reasonable event 
energy reconstruction only for quasi-elastic (QE) neutrino interactions. QE interactions 
produce just one charged particle above the Cerenkov threshold, resulting in just one 
Cerenkov ring in the detector. The two-body kinematics allows a precise reconstruction 
of the incident neutrino energy from the measured momentum of the outgoing lepton 
and the known direction of the incoming neutrino. E.g., for a QE i/^ event one has 

1 (M2 -ml) + 2EMn - Ml 

£jTcc ^ _J__P f^' p. n n (27) 

" "2 -E^ + Mn+p^, cos 9^ ^ ' 
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where Mp, Mn, m^, E^^ p^, cos6'^ are the proton, neutron, muon masses, muon energy, 
momentum and angle with respect to the incoming neutrino direction. Single ring events 
produced by non-QE interactions could bias the energy reconstruction because in their 
case the two body kinematics does not hold. For this reason it is important to precisely 
know the non-QE/QE ratio in the neutrino interactions. 

Other detector technologies, such as liquid argon TPCs, could suppress detector 
backgrounds to a negligible level and efficiently reconstruct neutrino interaction events 
of any multiplicity, and are very promising for future neutrino oscillation experiments. 

A beam configuration used to optimise long baseline experiments with the goal of 
measuring 6'i3 is the off-axis configuration, a concept that was ffistly proposed in |157] . 
According to the two-body vr-decay kinematics, all the pions above a given momentum 
produce neutrinos of similar energy at a given angle ^ 7^ with respect to the direction 
of the parent pion (contrary to the 6' = case where the neutrino energy is proportional 
to the pion momentum). 

An off-axis configuration offers several advantages for a long baseline experiment 
optimised to detect sub- leading z/^ — oscillations. First, the off-axis neutrino flux at 
the desired energy is higher than in an on-axis conflguration. Second, the neutrino flux 
at higher energies than the oscillation maximum is greatly reduced, with a consequent 
reduction of backgrounds due to 7r° generated by neutral current events. Both of these 
effects are clearly visible in comparing the MINOS (on-axis) and NOz/A (off-axis) spectra 
in flg. IHl Third, the intrinsic background due to beam v^. is less around the oscillation 
maximum (about 0.4% in T2K) since are mostly generated by three body meson 
decays that have different kinematics than the two body pion decays. These advantages 
overcompensate the smaller total neutrino flux (even when multiplied by the oscillation 
probability) of the off-axis conflguration. 

In an appearance experiment the influence of systematic errors is very much reduced 
if compared to the reactor experiments. Nevertheless a close detector is necessary in 
order to precisely measure the flux of v^^ that are going to oscillate into in the 
far detector, to measure interaction cross sections of signal and background processes, 
and to measure the backgrounds for the signal Ve- As extensively discussed in jl58] 
the task of the close detector is not just a background subtraction, since for instance, 
factors as the QE/non-QE ratio can spoil the experimental sensitivity. As it will be 
described in the following the T2K experiment renounced to have a replica of the far 
detector in the close position, and designed a very sophisticated close detector capable 
of precisely measuring the exclusive neutrino cross sections that play a role in the signal 
and backgrounds in the far detector as well as the flux of v^. NOz/A instead is planning 
to have just a replica of the far detector in the close station, but it shouldn't be missed 
the fact that in the same beam line (on-axis) will soon operate experiments specialised 
in measuring neutrino cross sections, like MINERz/A |159] . 

4-2.1. T2K The T2K (Tokai-to-Kamioka) experiment [2^ will use a high intensity 
off-axis (2.5°) neutrino beam, with a peak energy of 700 MeV, generated by a 30 GeV 
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Figure 11. Integrated exposures for the T2K (top) and NOz^A (bottom) experiments. 
The T2K beam power plan is based on [156) . The NOj^A exposure has been normalised 
to nominal "NOz^A years" , which corresponds to 6 x 10^" protons-on-target (pot) and a 
15 kt detector. The shown curves take into account the plan for the available detector 
mass as a function of time, as well as the available pot |160j . 

proton beam at J- PARC (Japan Proton Accelerator Research Complex) fired to the 
SuperKamiokande detector, located 295 km from the proton beam target. The design 
intensity of the J- PARC proton beam is 3.3 ■ 10^"^ protons/pulse with a repetition rate 
of 3.3 s, corresponding to a beam power of about 750 kW. The experiment is expected 
to collect data for about 10'' s/year, performances are normalised to a 5 year run at the 
nominal power. Fig. [11] shows the integrated exposure as a function of time using the 
expected beam power shown in |156] . According to this plan, the nominal exposure will 
be reached around 2017. 

The schematic view of the T2K neutrino beam line is shown in fig. [12] (left). The 
primary proton beam interacts in a helium-cooled graphite target placed inside the 
first horn. The beam optics includes three horns running at 320 kA with a maximum 
magnetic field of 2.1 T. The decay volume is 94 m long after the target region with a 
variable cross section starting at 2.2 m (W) at 2.8 m (H) and increasing to 3.0 m at 4.6 m 
at the far end, filled with inert helium gas held at 1 atm. A beam dump, constructed 
of cooled graphite and copper blocks is placed at the end of the decay volume. 

A sophisticated near detector complex has been built at a distance of 280 m from 
the target. This complex has two detectors: one on-axis (INGRID) and the other 
off-axis (ND280). This off-axis detector (fig. [121 right) is a spectrometer built inside 
the magnet of the former experiments UAl and NOMAD, operating with a magnetic 
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Figure 12. Left: the layout of the T2K beam hne, showing the location of primary 
proton beam line, target station, decay volume, beam dump, muon monitors and near 
neutrino detectors, [161 . Right: sketch of the T2K ND280 off-axis near detector. 



field of 0.2 T. It includes a Pi-Zero detector (POD), a tracking detector made by three 
time projection chambers (TPCs) and two fine grained scintillator detectors (FGDs) 
acting as an active target, a 4n electromagnetic calorimeter (Ecal), and a side muon 
range detector (SMRD). To reduce systematic errors both the POD and the FGDs 
incorporate water targets (the same material of the far detector). Neutrino rates in 
the close detector will be about 160000 z/^ (3200 z/g) interactions/ton/yr at the nominal 
beam intensity of 0.75 MW-10^ s. 

ND280 is expected to calibrate the absolute energy scale of the neutrino spectrum 
with 2% precision, measure the non-QE/QE ratio at the 5-10% and monitor the neutrino 
flux with better than 5% accuracy. The momentum resolution of muons from the 
charged current quasi-elastic interactions (CCQE) should be better than 10%. The 
Uf, fraction should be measured with an uncertainty better than 10%. A measurement 
of the neutrino beam direction, with a precision better than 1 mrad, is required from 
the on-axis detector. 

The close detector location, 280 m from the target, is too close to the decay tunnel 
to have a neutrino flux identical to the far detector flux. Indeed differences as big as 
50% are expected between the two fluxes, reducing the capability of the close detector 
of reducing systematic errors independently from any Monte Carlo simulation. A more 
distant close detector station could attenuate the near-far neutrino flux differences, but 
for the moment is not foreseen by the experiment. 

A fundamental tool to better control the neutrino flux is the ongoing measurement 
of the hadro-production in a T2K target replica by 30 GeV protons by the NA61 
experiment at CERN |162] . Following the succesful example of the measurements done 
by the HARP experiment [163] for the K2K [161] and for the MiniBooNE [165] targets, 
NA61 data should greatly improve the accuratness of the neutrino beam line simulation. 

Signal events are detected by the SuperKamiokande detector with an efficiency of 
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Near 
Detector 

Figure 13. Sketch of the NOz/A detectors. The insert figure shows the plane 
arrangement of each detector. From [166] . 

45% and with a detector background contamination smaller than the intrinsic beam Ue 
contamination (roughly in the 2:3 ratio, see also fig. [TOj) . 

4-2.2. NOvA The NOi/A experiment [23] will run at an upgraded NuMI neutrino 
beam expected to deliver 6.5- 10^*^ pot/year, corresponding to a beam power of 700 kW, 
generating a neutrino beam with an average energy ~ 2 GeV and a Vf. contamination 
less than 0.5%. 

The far detector, placed at baseline of 810 km, 14 mrad (0.8°) off-axis, will be a 
"totally active" tracking liquid scintillator, constructed from liquid scintillator contained 
inside extruded PVC cells. Each cell is 3.9 cm wide by 6.0 cm deep and is 15.5 meters 
long, they will be arranged in an alternating plane structure composed of vertical 
and horizontal cells. Scintillation light will be guided to APD photo-detectors using 
wavelength shifting fiber. NOz/A will be located near Ash River and will have a total 
mass of 15 kilotons (385000 cells) and be 15.7 meters wide, 15.7 meters tall, 78 meters 
long (the only non-cylindrical detector in this generation of experiments). It will run on 
surface, with a modest overburden of 3 m of concrete. It is expected to detect z/g signals 
with an efficiency of 26% (including acceptance) and to reduce the detector background 
rate to a level comparable to the rate from the intrinsic beam z/e contamination. The 
close detector will be a 215 ton replica of the far detector, placed 14 mrad off the NuMI 
beam axis at a distance of 1 km from the target. A sketch of the NOz/A detectors is 
shown in fig. [131 

NOz/A plans to alternate between neutrino mode and antineutrino mode: the 
focus of the experiment is to provide data on the neutrino mass hierarchy, where 
NOz/A has a clear advantage with respect to T2K thanks to the longer baseline, see 
e.g. [167] for a recent sensitivity study, and [168til70j for discussions of the T2K/N0z/A 
complementarity/synergy. A possible scenario [160] on the time evolution of neutrino 
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and antineutrino exposure is shown in fig. [TTl which takes into account the build-up of 
the detector mass as well as the available beam power. As a second phase, the NuMI 
beam intensity could be increased to 1.2 MW ("SNuMI") or to 2.3 MW ("Project X") 
in case the new proton driver of 8 GeV/c and 2 MW will be built at FNAL. 



5. Phenomenology and sensitivity estimates for upcoming experiments 
5.1. Complementarity of reactor and superbeam experiments 

A natural classification of the 5 upcoming experiments is the distinction between reactor 
and accelerator (superbeam) experiments. Apart from the vastly different experimental 
configurations and challenges (see section HI) there is also an important difference from 
the point of view of oscillation physics. While reactor experiments search for the 
disappearance of electron antineutrinos, accelerator experiments look for the appearance 
of electron (anti) neutrinos in a beam initially composed mainly of muon (anti)neutrinos. 

Let us now discuss the relevant oscillation probabilities. We denote with L,E^,V 
the distance from the neutrino source and the detector, the neutrino energy, and the 
matter potential, respectively, with V = ^J^GpNe where iVe is the electron density along 
the neutrino path. Furthermore we introduce the abbreviations 

Am^L , 1EV Am^ ^ ^ 

AE Am|i Am|i ^ ' 

Note that the signs of all three of these quantities depend on the neutrino mass hierarchy 
encoded in the sign of Am|]^. 

In the case of reactor experiments with L < 2 km the matter effect is negligible 
and the survival probability is given by 

Pee ~ 1 - sin^ 2^13 sin^ A - a^A^ sin^ 2^i2 . (29) 

As mentioned in the discussion of the CHOOZ experiment in section 13.11 the last 
term due to Am^i is suppressed by the small solar mass squared difference and can 
be neglected as long as sin^ 26^13 ~ 10"'^, and in this case Pee is independent of the 

solar parameters 612 and Amgi. Furthermore, the survival probability eq. |29]does neither 
depend on ^23, on the CP phase 5, nor on the sign of Amg^. Hence reactor experiments 
provide a "clean" measurement of ^13, free of correlations with other parameters, apart 
from |Am|^|, which, however, is relatively precisely known |15Hll52j . 

Let us now move to long-baseline appearance experiments. For typical neutrino 
energies in the GeV range and baselines larger than 100 km one cannot neglect 
the matter effect. However, for baselines smaller than a few 1000 km it is a good 
approximation to assume a constant matter density, given by the average density along 
the neutrino path. In that case a rather useful expression for the P^^^^y^ = P^e 
appearance probability can be obtained by considering terms up to second order in 
the small quantities S13 and a |171H173j : 

sin^(A - 1)A 

P^e^sm^2gi3 4 
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Figure 14. Exemplary fit results for Double Chooz, T2K, NOi/A, Daya Bay, and 
the combination. Shown are fits in the 613-6 plane assuming sin^ 2^13 — 0.1 and 
d — 7r/2 (upper row) and S — 37r/2 (lower row). A normal simulated hierarchy is 
assumed. The contours refer to Icr, 2cr, and 3<t (2 dof). The fit contours for the right 
fit hierarchy are shaded (coloured), the ones for the wrong fit hierarchy are shown as 
curves. The best-fit values are marked by diamonds and boxes for the right and wrong 
hierarchy, respectively, where the minimum for the wrong hierarchy is explicitly 
shown. Reprinted from Ref. \IQ71, Copyright (2009), with permission from JHEP. 



.X sin sin(A - 1)A 

+ a sin 26113 sin 26^12 sin 26^23 cos( A + 6) — ^- 

yi yi — 1 

9 0^0 sin^ AA , , 

+ sin^2gi2 4 ^2 • (30) 

This equation holds for neutrinos; for anti-neutrinos change 6 — > —b and V — )■ —V. The 
term in the first hne of eq. [30] is similar to a two-flavour oscillation probability, apart 
from the s\^ factor which controls the fraction of the flavour participating in the 
Vf. — oscillations. This term dominates for large ^13. The term in the second line of 
eq. [30] corresponds to an interference term between oscillations with Am|]^ and Am|^. 
It depends on the CP phase 8 and is responsible for CP violation. The term in the last 
line is independent of 6'i3 and describes oscillations with Am^i. It can be neglected for 
sin2 2^13 > 0.01. 

Obviously, the parameter dependence of the appearance probability is much more 
complicated than the one for disappearance in reactor experiments, as it depends 
on all 6 oscillation parameters. If information on ^13 is to be extracted from an 
appearance measurement the correlations with the parameters 5 and sgn(Am3]^) are 
especially important |174j . The difference between reactor and super beam experiments 
is illustrated in fig. [TH taken from |167j . It shows how typical fits in the ^i3-(5 plane 
would look like if ^13 was large (sin^ 2^13 = 0.1) and b was close to maximal CP violation 
5 = (upper rows) and 5 = 37r/2 (lower rows). Here NH has been assumed to generate 
the "data". Using this "true" hierarchy in the fit, the coloured regions are obtained. 
When the data is fitted with the "wrong" hierarchy, i.e., with IH in this case, the regions 
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Figure 15. Discovery potential of the five upcoming experiments in the plane of 
sin^ 2013 and S expected in 2018, see section 15.31 for our assumptions on exposure. 
To the right of the curves a non-zero value of 6*13 can be established at 3a. For the 
beam experiments we show normal (solid) and inverted (dashed) hierarchies, while 
reactor experiments are independent of the hierarchy. The four lines for Daya Bay 
correspond to different assumptions on the achieved systematic uncertainty, from 
weakest to strongest sensitivity: 0.6% correlated among detector modules at one site, 
0.38% correlated, 0.38% uncorrelated among modules, 0.18% uncorrelated. 



delimited by the curves are obtained. 

The figures show the characteristics of the different classes of experiments: The 
reactor experiments do not depend on 6, and the wrong fit hierarchy coincides with 
the right hierarchy. For T2K, which is simulated with neutrino running only, there 
is some dependence on 6, but the correlation between 6 and ^13 cannot be resolved. 
We observe the typical S-shape of the allowed regions, emerging from the cos (A + 6) 
term in eq. |30l The wrong hierarchy contours are slightly shifted, but the minimum 
is close to zero. NOz/A, on the other hand, has both neutrino and anti-neutrino 
running in the simulation, which means that the correlation can, at least in principle, 
be resolved. The wrong hierarchy can in some cases be excluded because of larger 
matter effects. In the combination of the experiments, the combination between Daya 
Bay and the beams allows for a substantial reduction of the allowed parameter space 
due to almost orthogonal measurements |151lll52lll75] . In the most optimistic cases, 
the mass hierarchy can be determined at 3a confidence level, and maximal CP violation 
can be established at relatively modest confidence as well. However, note that these 
optimistic cases represent only a very small fraction of the parameter space, see |167] 
for a discussion of the possibilities to measure CP violation and the mass hierarchy with 
these experiments. 

Fig. [15] shows the di^ discovery reach of the five upcoming experiments expected 
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Figure 16. Illustration of the luminosity scaling of the sin^ 26'i3 sensitivity at 90% CL 
of a fictitious reactor experiment. The horizontal axis shows the number of events in 
the far detector at a distance of 1.5 km from the reactor. We show the sensitivity 
for statistical errors only (red), as well as various values for the relative normalisation 
uncertainty between near and far detectors (blue curves with labels). The magenta 
and green curves show the sensitivity obtained for using only shape and only rate 
information, respectively. 

in 2018. Details about our exposure assumptions are given in section 15.31 It is clear 
from the figure that the discovery potential of the appearance experiments strongly 
depends on the CP-phase as well as on the neutrino mass hierarchy. We observe that 
the inverted hierarchy gives a weaker sensitivity. Hence, in case no appearance signal is 
found the final ^13 limit will be set by the IH. Note that this is completely analogous to 
the present information on ^13 from the MINOS z/g appearance data, compare figs. [15] 
and |3] (left). The fact that the IH has a weaker sensitivity is a consequence of the 
matter effect. It can be understood by considering the {A — 1) terms in eq. [301 Since 
MINOS as well as T2K operate only with neutrinos we have > 0, and A is positive 
(negative) for NH (IH). Hence, for NH (IH) the factors {A — 1) in the denominator 
lead to enhancement (suppression) of the oscillation probability, and therefore IH gives 
a weaker limit. The different shape of the IH curve for NOz/A results from the anti- 
neutrino running included in the NOi/A run plan. As evident from the figure, reactor 
experiments are neither sensitive to the value of 6 nor to the mass hierarchy. 

5.2. On statistical errors and systematical uncertainties 

Apart from the differences in the oscillation physics between reactor and super beam 
experiments there is also an important difference in how statistical and systematic errors 
influence the final sensitivity. 

Fig. [16] shows the sensitivity to sin^ 26'i3 as a function of the number of events for 
various assumptions on systematical uncertainties for a fictitious reactor experiment 
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with one reactor core and two detectors, at 0.15 and 1.5 km. The statistics-only 
sensitivity displays the expected 1/\/N scaling. The blue curves correspond to 
different values of the relative normalisation uncertainty between the two detectors. 
Turning on this uncertainty deteriorates the sensitivity above a characteristic exposure. 
Remarkably, at very high exposures again a 1/ a/ZV scaling is recovered |152] . This 
is a consequence of the information coming from the shape of the spectrum. The 
magenta line in the figure shows the limit from only shape information, with leaving 
the total number of events unconstrained in the fit. In reality also the spectral shape 
will suffer from some systematical uncertainty, which at some point will cut off again 
the 1/\/N scaling, see e.g., Refs. |1381I152] for a discussion. If only the total rate 
without any spectral information is used, the final sensitivity is just given by the over- 
all normalisation uncertainty, as illustrated by the green curve in fig. [16l To relate the 
three reactor experiments Double Chooz, RENO, and Daya Bay with fig. [16] we recall 
the expected events per year given in tab. [2] as (1.5, 3, 10) x 10^ respectively. Note, 
however, that the figure is for a fixed far detector baseline of 1.5 km and one single 
reactor, which does not correspond to any of the three experiments, and therefore the 
figure cannot be applied exactly to the specific experimental configurations. 

We note that especially the high-statistics experiment Daya Bay, with let's say 
3 X 10^ events after 3 years, is quite sensitive to the achieved systematic. This is 
illustrated in fig. [151 where we show the Daya Bay sensitivity for different choices 
on systematics. The four lines correspond to various assumptions on the relative 
normalisation uncertainty of the 8 detector modules. For the most conservative limit 
we assume the same uncertainty as claimed by Double Chooz, 0.6%, and take this error 
correlated between detector modules at each detector site. The Daya Bay "baseline" 
systematics is 0.38%. The two dashed curves correspond to this value assuming it either 
correlated or uncorrelated among detector modules at one site. The most aggressive 
curve shown in fig. [TSl assumes the "goal" value of the systematics of 0.18%, uncorrelated 
between all detector modules. Let us stress that of course systematics are also crucial 
for Double Chooz and RENO. In those cases we have adopted a systematical error of 
0.6%, as stated in the proposals. The final Double Chooz and RENO sensitivities will 
depend in a similar way on the assumed systematics value as exemplary shown for the 
case of Daya Bay in the plot. 

Let us now discuss the impact of systematics and statistics on the 6*13 measurement 
in super beam experiments. As an example we consider a sightly modified version of 
the T2K experiment, where, e.g., both neutrino and anti-neutrino running is assumed 
with a ratio of 1:3. A detailed description can be found in |158j . Fig. [17] (left) 
shows the smallest value of sin^ 2^13 which can be distinguished from Oi^ = as a 
function of the luminosity, assuming two representative values for the CP phase which 
correspond roughly to the best and worst sensitivity. The first observation is that for an 
exposure corresponding roughly to the nominal T2K exposure (marked by the vertical 
line labeled "T2K") systematics have only a small impact, since this measurement is 
largely dominated by statistics. Numerically, the sensitivity at a T2K exposure decreases 
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Figure 17. Left: Sensitivity to a non-zero 613 for a T2K-like experiment as a 
function of exposure for 5 — t:/2 and tt for a default choice of systematical errors 
and for statistical errors only (curves delimiting the shaded region). Furthermore, 
the sensitivity obtained without uncertainty on the intrinsic beam background (blue) 
and without an uncertainty on the NC background in the far detector (red-dashed) 
is shown. Right: Allowed region in the plane of sin^ 26*13 and 6 for two example 
choices for the input values marked by stars in the figure. Allowed regions are shown 
for all combinations of statistical errors only, systematics included, all other oscillation 
parameters fixed, and free (where for the solar parameters present errors are assumed). 
The sign(ATO|]^) degeneracy is neglected, and ^2™° = 7''/4. Reprinted from Ref. }158j . 
Copyright (2008), with permission from JHEP. 

from sin^ 2^13 = 0.0167 to 0.0172 for 6 = 7r/2, and from sin^ 2^13 = 0.0206 to 0.0214 for 
5 = n. 

Only for much larger exposures systematics have a non-negligible impact on the 
9i3 discovery reach. According to the discussion related to eq. [261 the uncertainty on 
the background is the most relevant systematics. Its impact is controlled by the ability 
of the near detector to predict the background in the far detector. Fig. [T7] shows also 
curves assuming a perfectly known i/g beam background, and no uncertainty at all on 
the background {i.e., fixing the beam contamination as well as the NC background). 
If the total background is fixed the sensitivity is close to the pure statistics case. It is 
interesting to note that for the two examples of 6 shown in the figure the importance of 
beam and NC backgrounds is different. This is an effect of the spectral shapes of the 
signal relative to the background, since the spectrum of the signal depends on the value 
of 6, and also beam and NC backgrounds have rather different shapes. 

Fig. [T7] (right) shows the allowed region in the space of sin^ 2^13 and 6 obtained by 
the same modified T2K-like configurationlJi] The impact of systematics is small, though 
not negligible in this case. Furthermore, we observe that the uncertainty on the other 
oscillation parameters has a sizable impact on the allowed region. This effect comes 
entirely from the atmospheric parameters Am^^ and 623. Apparently the disappearance 

II Note the different shapes of the regions from fig. [IT] (right) and the corresponding panels of fig. [HI 
The reason is the anti-neutrino running assumed for fig. 1171 
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channel does not provide enough accuracy on these parameters to avoid an effect on the 
6'i3 determination. For the solar parameters the accuracy from present data is sufficient 
to eliminate any effect on the results shown in the figure. 

5. 3. Sensitivity predictions for the upcoming experiments 

We now discuss the sensitivity of the different experiments to ^13, using two different 
performance indicators: the 613 sensitivity limit and the 613 discovery potential. The 
^13 sensitivity limit describes the ability of an experiment to constrain ^13 if no signal is 
seen. It is basically determined by the worst case parameter combination which may fake 
the simulated ^13 = 0. The sensitivity limit does not depend on the simulated hierarchy 
and 6, as ^^13 = is used to calculate the fake data. For a more detailed discussion, see 
appendix C of [138] . The 6*13 discovery potential is given by the smallest true value of 
6^13 > which cannot be fitted with 6'i3 = at a given CL. Since the simulated 6'i3, 6, 
and hierarchy determine the simulated rates, the 6*13 discovery potential will depend on 
the values of all these parameters chosen by nature. On the other hand, correlations 
and degeneracies are of minor importance because for the fit 613 = is used. 

For reactor experiments both measures (sensitivity limit and discovery potential) 
are very similar, since statistics are Gaussian and the oscillation physics is simple. For 
beam experiments the smallest ^13 discovery potential for all values of 5 and mass 
hierarchies (risk-minimised ^13 discovery potential) is often similar to the ^13 sensitivity 
limit. However, notable deviations from this rule occur due to Poisson statistics as 
well as more complicated oscillation physics implying correlations and degeneracies. In 
particular, we find that for low exposures the worst case discovery potential of the beams 
is somewhat better than the limit in case of no signal at the same confidence level. 

We are going to show the evolution of the 6*13 sensitivity with time, assuming the 
following schedules for the experiments based on up-to-date information. For Double 
Chooz the far detector starts data taking in June 2010, while near detector data will 
become available beginning of 2012. RENO starts data taking with both close and 
far detectors end of 2010. For Daya Bay we do not consider the periods where only 
the near detectors are available in 2011, but consider the start of the full experiment 
with all detectors at the end of 2011. For all reactor experiments we assume that all 
reactors are at nominal power all time according to the event numbers per year given in 
table O For T2K and NOi/A we use the exposures as a function of time shown in fig. HH 
assuming that protons on target are uniformly distributed along the year, not taking 
into consideration the specific schedules of the accelerators. In all cases we assume that 
results are available instantaneously with data. 

The simulations of T2K and NOi^A are performed with the GLoBES software [133[ 
I134j based on the experiment definitions developed in |167j and available at the GLoBES 
web-page. Modifications to the T2K simulation have been introduced due to recent 
updates on efficiencies (it now reproduces fig. [TOl) . For the reactor experiments an 
independent code has been developed, allowing for an arbitrary number of reactors and 
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Figure 18. Evolution of tlie ^13 sensitivity limit as a function of time (90% CL), i.e., 
the 90% CL limit which will be obtained if the true ^13 is zero. The four curves for 
Daya Bay correspond to different assumptions on the achieved systematic uncertainty, 
from weakest to strongest sensitivity: 0.6% correlated among detector modules at one 
site, 0.38% correlated, 0.38% uncorrelated among modules, 0.18% uncorrelated. 

detectors. Various systematics are included with proper correlations between detectors 
and reactors, as well as backgrounds from accidental, fast neutrons, and cosmogentics 
according to the numbers provided in the respective proposals. Reactor fluxes and their 
uncertainties are included following |176j . 

The sensitivity limit time evolution is shown in fig. [181 We observe that 
the global sensitivity limit will be dominated by reactor experiments. If the assumed 
schedules for the reactor experiments are achieved. Double Chooz and RENO will each 
dominate the limit between 0.5 and one year. As soon as Daya Bay becomes operational 
it will have the best limit (as well as discovery potential) among the reactor experiments, 
thanks to the large exposure. After 5 years of running the limits at 90% CL will be 
sin^26'i3 = 0.026 and 0.019 for Double Chooz and RENO, respectively, whereas the 
Daya Bay limit will range from 0.012 to 0.006, depending on the systematics. In fig. [18] 
we show the Daya Bay limit under various assumptions about the relative normalisation 
uncertainty of their detectors, as discussed in section 15.21 Even for the most pessimistic 
assumption, Daya Bay will set the final limit on 613 in the case no finite value will be 
discovered. 

In case of no signal, the 61^ limit from beam experiments suffers from the 
marginalization over the CP phase and the mass hierarchy, as discussed in section 15.11 
This situation is very different in case of the discovery potential, since there a favourable 
value of 6 can greatly enhance the sensitivity of the appearance experiments. The 613 
discovery potentials are shown in fig.[T9]as a function of time. For the beam experiments, 
the dependence on the true value of 6 is reflected by the interval between the solid curves 
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Figure 19. Evolution of the 6*13 discovery potential as a function of time (3ct CL), i.e., 
the smallest value of ^13 which can be distinguished from zero at 3cr. We assume the 
normal and inverted simulated hierarchies in the top and bottom panels, respectively. 
The bands for the beams reflect the (unknown) true value of S. For the dashed curves 
S — has been fixed. The four curves for Daya Bay correspond to different assumptions 
on the achieved systematic uncertainty, from weakest to strongest sensitivity: 0.6% 
correlated among detector modules at one site, 0.38% correlated, 0.38% uncorrelated 
among modules, 0.18% uncorrelated. 



for a given time (shaded regions). The dashed curves for T2K and NOz/A correspond to 
a fixed value for the CP phase of 5 = 0. The reactor experiments are not affected by 
the true 6; the various curves for Daya Bay again correspond to the different assumptions 
concerning systematics as described above. There is a small dependence on the true mass 
hierarchy for the beam experiments, compare top and bottom panels, where for IH the 
sensitivity is slightly worse because of the suppression of the oscillation probability for 

% Evolution of sensitivities under this condition have been shown recently in |1 7711178] . 
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Figure 20. Evolution of the ^13 discovery potential as a function of time (Scr CL) for 
NH, showing the global sensitivity reach. The bands for the beams and the global reach 
reflect the (unknown) true value of S. For Daya Bay we have assumed a systematical 
uncertainty of 0.38% correlated among detector modules at one site. 



neutrinos by the matter effect. 

The comparison of Figs. ITQ] and fTSl shows that suitable values of S may significantly 
improve the discovery potential of beams compared to their sensitivity limit. Indeed, 
T2K may discover ^13 for smaller 613 than Daya Bay in a significant fraction of the 
parameter space, depending on the achieved systematics in Daya Bay. The NOz/A band 
becomes more narrow due to the complementary information from the anti-neutrino 
running, with the clear disadvantage of being somewhat late. 

In figure [20] we illustrate how the world sensitivity to 613 could look like under the 
assumptions of the above schedules and that at each point in time a combined analysis 
of all available data is performed. The discovery reach will be set roughly by the optimal 
sensitivity of T2K, where the reactor experiments play an important role in providing 
sensitivity for the values of S unfavourable for T2K. This plot nicely illustrates the 
interplay between reactor and beam experiments and shows that the global reach can 
be enhanced significantly if experiments of both types are available simultaneously with 
comparable sensitivities. 

Let us comment also on the dependence of the sensitivities on the assumption 
sin^ 023 = 0.5. This is relevant only for beam experiments, since the survival probability 
in the reactor experiments is independent of 623, see eq. [291 In the range 0.4 < sin^ 6*23 < 
0.6 (corresponding roughly to the current 2a allowed range [H]) the effect for the beam 
experiments is small |167] . The sensitivity to ^13 becomes somewhat better the larger 
sin^ 023, which follows from the first term in eq. [301 In these cases the choice sin^ 623 = 0.5 
corresponds to the "average" situation. 

Note that this discussion is based on the unitary standard three-flavour oscillation 
framework. If the search for new physics is taken into account, different reactor 
experiments, or reactor experiments and superbeams, may imply different information 
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and therefore be very complementary; see, e.g., Refs. |179[I180] . 

5.4- On the interpretation of sensitivities for future experiments 

A widely used procedure to calculate sensitivities for future experiments is to assume 
some input values ("true" or "simulated" values) for the oscillation parameters for 
which the predictions for the observables in a given experiment are calculated without 
statistical fluctuations. Then these predictions are used as "data" and a statistical 
analysis of these data is performed to see how well the input values for the parameters can 
be reconstructed by the experiment. This procedure, denoted by "standard" procedure 
in the following, should give the sensitivity of an "average" experiment, where "average" 
lacks a precise definition. The results of the previous subsections are obtained by this 
method, and most of the sensitivity studies for future neutrino oscillation experiments 
in the literature are based on this procedure. In particular the GLoBES long-baseline 
experiment simulation software |133lll34j is designed for this method. 

Ref. [181] clarifies the correct interpretation of such sensitivities in the context 
of oscillation experiments. Monte Carlo simulations of the Double Chooz and T2K 
experiments are performed in order to address the following question in a frequentist 
framework: Given a true value 0/6*13, what is the probability that the hypothesis 6*13 = 
can be excluded at a certain confidence level? This generalises the usual sensitivity limits 
to a well defined statistical statement and allows also a precise definition of the "average 
experiment". For example one may define the sensitivity of an average experiment as 
the value of ^^3™ for which 613 = can be excluded with a probability of 50%. 

A large number of artificial data sets is generated to calculate the actual distribution 
of the statistics used to decide whether ^13 = should be rejected at a given confidence 
level. This allows to answer the question stated above within a well defined frequentist 
framework. Moreover one does not rely on questionable assumptions necessary in 
the standard procedure, for example issues related to the non-linear character of the 
parameters, the periodicity of the CP phase 6, the physical boundary sin^ 26*13 > 0, 
assuming standard x^-distributions, and the question of how many degrees of freedom 
to use for them. Note that the specific experimental configurations and parameters 
used for the simulations in |181j for Double Chooz and T2K differ slightly from the 
most up-to-date versions. 

The probability -Pdisc(«, 6*13) that ^13 = can be excluded at the 100(1 — a)% CL 
is given by 



Here A^q = x^(6'i3 = 0) — xLin is the difference in between the best fit point and 
^13 = 0, \{a) is the Ax^ value corresponding to the 100(1 — a)% confidence level, 
and /0^3(Axq) is the distribution of Axl for a fixed true value of 6*13. The distribution 
fe-^^i^Axo) and the values of A(a) have been obtained by Monte Carlo simulation. 




(31) 
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Figure 21. The probability to exclude the hypothesis 613 = at the 99.73% CL 
for T2K and Double Chooz as a function of the true value of sin^ 29i3. The two 
curves for T2K correspond to the true values 6 = 108° and 288°. The vertical lines 
show the corresponding "standard" sensitivities. The dashed curves correspond to 
the probability Pdisc calculated in the Gaussian approximation according to eq. |32l 
Reprinted from Ref. |181) . Copyright (2007), with permission from Elsevier. 



Fig. [2T] shows the probabihty Pdisc to exclude the hypothesis 6^13 = at the 
99.73% CL for T2K and Double Chooz as a function of the true value of sin^2^i3. For 
each true value 3 x 10^ data sets have been simulated. For T2K the two values chosen 
for S correspond roughly to the best and worst sensitivity. The vertical lines in the plot 
show the standard sensitivities calculated from the condition Ax^ > 9 without statistical 
fluctuations. One observes that for Double Chooz the standard sensitivity corresponds 
indeed with good accuracy to Pdisc = 50%, as expected for an "average" experiment. For 
T2K the discovery probabilities corresponding to the standard sensitivities are actually 
slightly higher, around 60%. 

The dashed curves shown in fig. |2T]are obtained assuming a Gaussian measurement 
of sin^26'i3. In this case Pdisc can be obtained in terms of the error function in the 
following way. Assuming that a; is a Gaussian variable with standard deviation a the 
hypothesis x = can be excluded at the 99.73% CL if the observed value x°^^ is bigger 
than 3a. On the other hand, the probability for x°^^ > 3cr as a function of the true 
value x*™^ is easily calculated as 

1 - ert I 

2 



^00 -1 

P [x"'"' > 3a I x'''"'] = I dx G{x; a) 



V2, 



a 



(32) 



'3a 

where G{x; x*'''^*^, a) denotes the normal distribution with mean x*''"^ and standard 
deviation a. 

The dashed curves in fig. |2T] have been obtained from eq. [321 by identifying 
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Probability to observe a non-zero at 99.73% CL in T2K 




Figure 22. Contours of the probability Pdisc to establish a non-zero value of O^i at 
the 99.73% CL for T2K in the sin^ 26''i5''=-(5*™'= plane. The dashed curve corresponds 
to the "standard sensitivity limit". Reprinted from Ref. |181) . Copyright (2007), with 
permission from Elsevier. 



sin^20i3 = X and by using for cr one third of the 99.73% CL sensitivity hmit from 
the standard procedure. One observes that for Double Chooz this approximation is 
excellent. Hence, in this case sin^ 26*13 can be considered indeed as a Gaussian variable 
and the probability Pdisc can be calculated from the standard sensitivity limit and 
eq. |32] without the need of a MC simulation. In contrast, for T2K some deviations 
from Gaussianity are visible (especially for 5*''"*^ = 288°). This is not unexpected, since 
in this case event numbers are small, background fluctuations are important, and the 
dependence of the observables on the parameters is much more complicated than in the 
case of Double Chooz. 

Contours of the probability Pdisc for the T2K experiment in the plane of sin^ 29^^"^ 
and 6^"^^^ are shown in fig. [221 -Pdisc has been calculated for a grid of 41 x 41 values 
and at each point in the grid 10^ data sets have been generated, leading in total to 
nearly 1.7 x 10^ performed fits. This figure is the generalisation of the usual sensitivity 
limit (shown as dashed curve) and for each true value of the parameters one can infer 
the probability that T2K can establish a non-zero value of ^13 at the 99.73% CL. As 
indicated already in fig. [2T] one finds that the standard sensitivity curve corresponds 
roughly to a discovery probability of 60%. The region where 6*13 > can be established 
with high probability, let's say greater than 99%, is found for sin^ 29f^^'^ > 0.0166 — 0.041, 
depending on the true value of 6. It is shifted with respect to the standard sensitivity 
limit to values of sin^ 2^13 larger by roughly a factor of 2. 

In summary, for Double Chooz (and reactor experiments in general) the Gaussian 
approximation is very well justified. The usually calculated sensitivity corresponds to 
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the performance of an average experiment (the discovery will be made with a probability 
of 50%), and the actual discovery probability can be estimated by a simple formula in 
terms of the error function. In the case of appearance beam experiments some deviations 
from Gaussianity are found. For T2K the standard sensitivity limits correspond to a 
discovery probability of about 60%. The results of |181] confirm that standard sensitivity 
limits provide a reasonable approximation for an average experiment, corresponding to 
a discovery probability of order 50%. However, one has to be aware of the correct 
interpretation of such limits. In general the region where a discovery can be made with 
high probability is significantly smaller than the one corresponding to the standard 
sensitivity limits. 

6. A subsequent generation of experiments 

If the upcoming generation of experiments discussed in this review cannot establish a 
finite value of ^13, experiments able to explore the region sin^ 2^13 < 10^^ will be needed. 
On the other hand, if a positive signal can be established this generation of experiments 
will not be sensitive enough to make a firm discovery (3cr or better) of leptonic CP 
violation and the neutrino mass hierarchy, see Fig. [14] and the extensive discussion in 
Ref. |167] . In this section we briefely outline some possibilities for such high precission 
oscillation facilities, more complete discussions can be found in [29||33]. 

A possible way to attack the searches of leptonic CP violation is to push 
conventional neutrino beams to their ultimate performances and build huge detectors, 
one order of magnitude, at least, bigger than the present detectors. As an example T2K 
could be upgraded, as already dehneated in the initial Lol [21], by pushing the power of 
the J-PARC Main Ring to 4 MW (2 MW seems nowadays more realistic) and building 
a water Cerenkov detector of about 500 kton fiducial volume, HyperKamiokande, more 
than 20 times bigger than SuperKamiokande, to be placed at the same distance and 
the same off-axis angle as SuperKamiokande. Subsequent developments of this project 
foresee the possibility of placing half of the detector at a longer distance, about 900 km 
in Korea, T2KK |182] . or to substitute it with a 100 kton liquid argon detector placed at 
658 km from J-PARC at a smaller off-axis angle of 0.5° |183j . The longer baselines would 
provide better sensitivities to the mass hierarchy and better control of degeneracies. In 
the United States a similar project designs a wide-band beam (WBB) |184] generated 
by the FNAL Main Injector and fired to a a 300 kton water Cerenkov detector placed at 
the DUSEL lab, 1290 km far away from FNAL (3^6 liquid argon modules of 20 kton 
are also taken in consideration). 

In Europe, superbeams have been proposed based on upgrades of the CNGS 
beam |1391I140] . or on a high power version of accelerators foreseen for a possible 
new injection chain of the LHC |185] as the SPL and the PS2. It appears very 
difficult |186] to set the performances of the CNGS beam to the level needed for 
searches of leptonic CP violation, that is 10 times the protons actually generated on 
target. The SPL superbeam |187j . based on the 3.5 GeV, 4 MW Superconducting 
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Proton Linac |188] . can generate a neutrino beam of about 300 MeV sent to a 500 kton 
water Cerenkov detector (MEMPHYS [HI]) placed in Frejus at 130 km from CERN. 
It would have excellent sensitivities for ^13 and leptonic CP violation |125] . while the 
mass hierarchy sensitivity would be limited by the relatively short baseline (the synergy 
with atmospheric neutrinos [73] would provide sensitivity to the mass hierarchy at large 
values of [125]). 

It has been proposed in |190] to generate a neutrino beam by a high power (1.6 MW) 
version of the PS2 accelerator, a 50 GeV synchrotron designed to run at 0.4 MW to 
serve as a component of the new injection scheme for the LHC Neutrinos could be then 
fired to a 100 kton liquid argon detector, placed at distances between 950 km or 1544 km 
or 2300 km. The distances correspond to the three underground labs of Sieroszowice in 
Poland, Slanic in Romania and Pyhasalmi in Finland, respectively, considered by the 
LACUNA FP7 Design Study [191j. As in the case of the WBB at DUSEL, this setup 
would measure neutrinos at the first and at the second oscillation maximum. Liquid 
argon is certainly the best candidate to fulfill the requirements of this configuration. 
This kind of configuration would have excellent performances in measuring the sign of 
Am|^ but a limited sensitivity for leptonic CP violation and the measurement of ^13. 

Sensitivities on 6*13 of these future superbeams could reach the sin^(26'i3) ~ 10^^ 
level, as illustrated in Fig. |23l Similarly also leptonic CP violation could be performed 
in the same range of ^13 values, see e.g., [2911301133] . 

The main limitations of conventional neutrino beams are the intrinsic 
contamination of the beam produced by the decay-in-flight of the kaons produced at the 
target and by the muons from the pion decays, and the uncertainty of the beam flux 
prediction due to the finite precision with which the hadron production cross sections 
can be known. The intrinsic limitations of conventional neutrino beams can be overcome 
if the neutrino parents are fully selected, collimated and accelerated to a given energy. 
This can be attempted within the muon lifetime (neutrino factory [25]) or within beta 
decaying ion lifetimes (beta beam [2^). In this way very pure, intense and precisely 
known neutrino beams can be designed. 

The beta beam [26] is a facility based on the decay in flight of /3-unstable ions, 
for reviews see [271 128]. They are ideal tools to study z/e — )■ transitions and their 
CP-conjugate without any intrinsic source of backgrounds. The original proposal of 
Ref. [26] was tuned to leverage at most the present facilities of CERN (the PS and 
the SPS) and it was based on ^He and ^^Ne as Ve and Ue sources, respectively. These 
ions are accelerated to 7 = 100 [195] and stored in a decay ring of about 7 km. With 
MEMPHYS as a far detector at Frejus, this set-up could outperform any superbeam 
configuration in terms of sensitivity to 6*13 and leptonic CP violation [125] . and these 
performances would be significantly enhanced if the SPL superbeam would be fired to 
the same far detector [125lll96j . Accelerating the same ions to 7 = 350 (an option that 
requires a 1 TeV accelerator) the performances of the beta beam would be significantly 
improved [19711198] . Various even more ambitious beta beam configurations are discussed 
in the literature, including very high 7 factors as well as exploring different isotopes. 
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Figure 23. 3a discovery sensitivity to sin 2^13 as a function of the fraction of all 
possible values of (Sep, computed for 10 years of data taking and 5% systematic errors 
(2% for the neutrino factory). The curve labelled WBB-WC refers to the FNAL- 
DUSEL setup [192j . the curves for the SPL superbeam, the 7 = 100 beta beam 
and their combination are taken from |125] . the curve for the 7 = 350 beta beam 
is taken from [193) . the low energy neutrino factory (LENF) from |194) . and IDS-NF 
corresponds to a hypothetical neutrino factory scenario with two 100 kton magnetized 
iron detectors at 4000 and 7500 km and an Emulsion Cloud Chamber detector to detect 
i>e Vt transitions at 4000 km taken from |33j . The plot has been adapted from a 
figure prepared by P. Huber for the EUROj^ WP6 study [33] . 



see [33] for a recent summary and references. 

Production, acceleration and stacking of high intensity muon beams for muon 
colhders have been envisaged thanks to their decays producing useful beams of and z/g 
(exploiting /i~ decays into e'Vf.i'^) or and z/g (/^^ decays into e'^Pe'^ii)- In a "Neutrino 
Factory" ^25j , muons are created from an intense pion source at low energies, their phase 
space is compressed to produce a bright muon beam, which is then accelerated to the 
desired energy and injected into a storage ring with long straight sections pointing in 
the desired direction. The neutrino factory design can be considered as strongly physics- 
motivated intermediate step towards a muon collider |199] . 

Assuming /i"*" in the decay ring, one looks for z/g — )■ z/^ oscillations due to the 
appearance of /x~ from CC events in the detector ("wrong sign muons"), which have 
to be separated from the bulk of /i+ ( "right sign muons" ) coming from unoscillated V^. A 
suitable detector to search for these transitions is a magnetized iron detector |200p201j . 
As pointed out in |202] , a neutrino factory would be an ideal tool to address CP violation 
in the leptonic sector, with outstanding performances compared with pion-based sources. 
The realization of the neutrino factory still represents a major accelerator challenge 



Ois-' phenomenology, present status and prospect 



46 



compared with neutrino superbeams. 

We report in Fig. [23] a comparison of sensitivities of some examples for the 
facihties described in this section. 

7. Conclusions 

In this review we have summarized the present status of the last unknown lepton mixing 
angle 9is and discussed the experimental progress to be expected within the next years. 
Although there some hints for a non-zero value of ^13 in current neutrino oscillation 
data, at present we can only quote an upper bound of sin^ 613 < 0.031 at 90% CL from 
global data. Best fit points obtained by different groups range from sin^^is ^ 0.01 to 
0.02, being consistent with 6*13 = at the 1 to 2a level, see Tab. [T]for a summary. 

Several new neutrino oscillation experiments will contribute significantly to our 
information on ^13 in the near future. The reactor experiments Daya Bay, Double Chooz, 
and RENO as well as the accelerator experiments NOz/A and T2K have been presented 
in section |H The time evolution for the sensitivity to ^^13 from these experiments is 
summarised in the figures [18] and [191 Within the next few years values of sin^ 26'i3 down 
to 10~^ will be probed, about one order of magnitude smaller than the current upper 
bound. These results will be of fundamental importance for the better understanding 
of the lepton sector, the problem of flavour and in particular for the question whether 
the values of the lepton mixing matrix indicate the presence of a flavour symmetry. 
The results of those experiments will be a corner stone for any subsequent neutrino 
oscillation facility aiming at the ultimate goals like the discovery of CP violation in the 
lepton sector or the determination of the neutrino mass hierarchy. 
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